Quantum criticality in a partially geometrically frustrated heavy-fermion system: 

CePdi_^Ni^Al 
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In the antiferromagnetic (AF) heavy-fermion system system CePdAl the magnetic Ce ions form a 
network of equilateral triangles in the (001) plane, similar to the kagome lattice, with one third of the 
Ce moments not participating in the long-range order. The Neel temperature Tn = 2.7K is reduced 
upon replacing Pd by Ni in CePdi_2:Nia,Al, with Tn — >■ for x — 0.144. At this concentration the 
specific heat C exhibits a C/T oc — log T dependence. Within the Ifertz-Millis-Moriya model of 
quantum criticality, this behavior might indicate 2D critical antiferromagnetic fluctuations, which 
might arise from the decoupling of the 3D magnetic order by frustration. On the other hand the 
simultaneous presence of Kondo effect and magnetic geometric frustration might entail a new route 
to quantum criticality. 



Quantum phase transitions (OPT) have received a lot 
of attention in recent years [l|-|3|. An issue of current 
interest is the role of dimensionality when a system is 
approaching a quantum critical point (QCP). A partic- 
ularly compelling case is given by heavy-fermion (HP) 
systems where the transition between a magnetically or- 
dered ground state and a paramagnetic Permi-liquid state 
is driven by the competition between intersite Ruderman- 
Kittel-Kasuya-Yosida (RKKY) and on-site Kondo in- 
teractions, which can be tuned by external pressure or 
chemical composition Unconventional quantum crit- 
ical behavior was observed in the canonical HP system 
CeCug-xAua; where quasi-2D magnetic fluctuations arise 
out of 3D long-range antiferromagnetic (AP) order, as 
shown by inelastic neutron scattering Q. Furthermore, 
the QCP is associated with unusual energy /temperature 
scaling of the dynamical susceptibility arising from crit- 
ical fluctuations Q . These findings prompted the de- 
velopment of alternative models of quantum criticality 
0, iSi] going beyond the Landau-Ginzburg- Wilson model 
of classical second-order phase transitions and its exten- 
sion to QPTs at zero temperature by Hertz P] , Millis [lo| 
and Moriya ll| (HMM model). In the model of local 
quantum criticality derived within the framework of ex- 
tended dynamical mean- field theory Q, 2D fluctuations 
are a prerequisite for the validity of the model. The is- 
sue of dimensional crossover in metallic quantum-critical 
magnets was discussed theoretically by Garst et al. [l^ . 
The issue of dimensionality was also addressed recently 
by comparing several systems of reduced dimensionality 
with cubic Ce3Pd2oSi6 at a magnetic QCP [l3| . 

Geometric frustration is a possible means to access 
quantum critical behavior in magnets on more general 
terms, with the possible appearance of spin liquids, e.g., 
in kagome lattices 
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Consequently, a more gen- 
phase diagram of quantum criticality was 
predicted |l5l4l7| . Hence it is of considerable interest 
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to study quantum criticality in HP systems with par- 
tial frustration. Pirst experiments were done recently 
on Ce2Pt2Pb, a realization of the Shastry-Sutherland 
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FIG. 1. (a) View on the ab plane of CePdAl. The Ce and 
Pd(2) atoms form a plane at z = 0, the Al and Pd(l) atoms 
form a plane at z = 1/2. (b) View on the a6-plane, show- 
ing one of the three symmetry-related magnetic structures. 
Only Ce-atoms are shown for clarity, (c) Kagome lattice for 
comparison. 



model, without however accessing quantum critical be- 
havior [l8| . 

In this paper we investigate the quantum critical be- 
havior of partially frustrated CePdAl, a HP system with 
a distorted kagome lattice, see Fig. [1] (a). In the hexag- 
onal ZrNiAl-type crystal structure of CePdAl, the mag- 
netic Ce ions form a network of equilateral corner-sharing 
triangles in the ab plane [3,[23l, see Fig.[l](b). The com- 
pound exhibits a strong magnetic anisotropy with the 
susceptibility ratio Xc/Xab ~ 14 which is attributed to 
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crystalline-electric- field and exchange anisotropics 2l|. 
CePdAl exhibits Kondo-lattice properties as evidenced 
by an increase of the electrical resistivity towards low 
temperature yielding a Kondo-temperature Tk « 5 K 
[23] . and orders antiferromagnetically below Tn = 2.7 K. 
A partially ordered magnetic state of CePdAl was first 
revealed by neutron diffraction measurements on pow- 
der samples [lO], where one third of the Ce moments 
do not participate in the long-range order. The mag- 
netic structure of CePdAl consists of three inequivalent 
Ce sites, with a magnetic ordering vector Q = (i r) 
with T w 0.35 being weakly temperature dependent [23j . 
Within a single kagome-like layer, one thus has ferro- 
magnetic chains coupled antiferromagnetically (see Fig.[T] 
(b)). A recent thorough ^"^Al NMR study indicated that 
the partially ordered state is indeed stable down to at 
least 30 mK (Ref.ial). 

The nearest-neighbor (nn) and next-nearest neighbor 
(nnn) Ce-Ce distances within the basal plane are re- 
sorted as dahi — 3.73 A and dab2 — 5.25 A respectively 
2o[ [25I . The interplane nn distance is d±^ = 4.24 A 
201 . Neglecting the interplane interactions the magnetic 
structure can be described in terms of a purely 2D model 
with nn ferromagnetic(Ji) and nnn AF ( J2) interactions 
within the ab plane [26[, see Fig. [1] (b). The competi- 
tion between Kondo and RKKY interactions was mod- 
eled in this approach by a T-dependent on-site parameter 
Ai (T). The effect of geometric frustration on this compe- 
tition was recently studied in detail by variational Monte 
Carlo simulations [i^l- Unfortunately, a model incorpo- 
rating the interplane coupling J±^ to be compared with 
the experimentally determined magnetic structure does 
not exist up to now. The QCP of CePdAl can be accessed 
by partial substitution of Pd by Ni [H HI]. We find 
that Tn of CePdAl is monotonically reduced to T = 
with increasing Ni content. At the critical concentration 
Xc ~ 0.144, the specific heat C as a function of temper- 
ature T varies as C/T cx log(ro/r), within the HMM 
model suggestive of 2D critical fluctuations. 

Polycrystalline samples of CePdi_a;Nia;Al were pre- 
pared by arc-melting appropriate amounts of the pure 
elements Ce(99.9) Pd(99.95), Ni(99.95), Al(99.999) 
under argon atmosphere with titanium gettering. To 
achieve homogeneity, the samples were turned over af- 
ter melting and solidifying and remelted several times. 
The total weight loss after preparation did not exceed 
0.5%. The CePdi-ajNij^Al samples were investigated in 
the as-cast state since annealing may cause a change 
of their structure [3l|. All samples were characterized 
by powder x-ray diffraction and atom absorption spec- 
troscopy (AAS). They revealed the single-phase ZrNiAl 
structure (P62m) of the parent compounds CePdAl and 
CeNiAl. AAS was used to determine the actual Ni 
concentrations x, that are quoted throughout this pa- 
per. The lattice constants a and c and the unit-cell vol- 
ume V approximately follow Vegard's law in the con- 




FIG. 2. (a) Specific heat C of CePdi-^Ni,;AI plotted as C/T 
versus T. (b) Specific heat C of CePdo.856Nio.i44Al plotted as 
C/T versus . (c) Magnetic entropy S versus T for x = Q. 



centration range investigated {x < 0.15). Specific-heat 
measurements were performed in a dilution refrigera- 
tor in the temperature range 0.05 < T < 2.5 K using 
the standard heat-pulse technique. A Physical Prop- 
erties Measurement System (PPMS, Quantum Design) 
was used to obtain data at higher temperatures for some 
samples. The dc magnetic susceptibility x "^^^ mea- 
sured at 0.1 T in the zero- field-cooled field- heated mode 
in a vibrating sample magnetometer (VSM, Oxford In- 
struments). A sample-dependent residual background 
contribution xo ~ 2.09- lO"'' ^s/f.u. T (x = 0) and 
2.09 • lO-^/zs/f.u.T {x = 0.144) independent of tem- 
perature was subtracted from the data. 

The specific-heat data are shown as C/T vs. logT in 
Fig. [21 The pure CePdAl compound (x — 0) exhibits 
a sharp anomaly at the Neel temperature Tn ~ 2.7 K 
in agreement with literature data 32|. The anomaly 
broadens and moves to lower temperatures with increas- 
ing Ni content indicating a suppression of the antifer- 
romagnetic (AF) transition. For the concentration x — 
0.144, the C/T vs. logT data follow a straight hue in 
Fig. [2] over almost two decades in temperature in the 
range 0.05 < T < 3 K, i.e., C/T = alog{To/T), with 
a = 0.705 J/mol and To = 11.7 K. The non-Fermi- 
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liquid behavior in the form of a logarithmic divergence 
of C/T versus T was previously observed in a number 
of systems, including CeCug-xAua; for x = 0.1 at the 
concentration-tuned QCP and samples with x > 0.1 at 
the pressure induced QCP 13)]. The HMM model [§-[11 
predicts for 2D antiferromagnets a logarithmic divergence 
of C/T towards low T near the QCP. Hence the C/T 
dependence of CePdi_a;Ni2:Al can in this model be ex- 
plained by 2D critical fluctuations, while the HMM pre- 
diction for 3D antiferromagnets, C/T oc 7 — aVT for 
T — 0, ist not compatible with our data, as shown in 
Fig [2] (b) . The magnetic entropy S is expected to be 
S — i?ln2 = 5.76J/molK for the twofold degenerate 
crystal-field groundstate of hexagonal CePdAl. Taking 
into account that one third of the Ce ions remain disor- 
dered down to lowest temperatures, one would expect an 
entropy S ~ |i?ln2 = 3.84J/molK to be recovered at 
T w Tn- Fig. [2] (c) shows that this is approximately the 
case for pure CePdAl. 

We use the specific- heat data C{T) plotted on a linear 
scale (not shown) defining Tn as the temperature where 
C(T) is a maximum for low x. For larger Ni concentra- 
tions, however, the small specific-heat anomaly associ- 
ated with magnetic order resides on a large T-dependent 
background Cbg, similar to the behavior of CeCue-xAux 
close to the QCP [s^]. Hence we assume that the spe- 
cific heat of the allegedly quantum-critical sample with 
X — 0.144 presents the 'background' that is subtracted 
from the data of the samples close to Xc and take Tn as 
the temperature of the maximum of the resulting AC{T) 
curve. Of course, we can not pin down the critical con- 
centration with absolute accuracy. If the critical con- 
centration would be slightly above 0.144, we estimate 
for X = 0.144 a maximum possible broad anomaly at 
Tn to ss 0.15 K for this sample. The resulting Tn(x) is 
shown in Fig. |3] (upper horizontal scale). Up to a; = 0.1, 
Tn(x) nicely follows a straight line as expected for the 
2D AF QCP in the HMM model. The deviations to- 
wards the QCP possibly signaling a restauration of three- 
dimensionality, will be discussed below. 

In order to compare the x dependence of Tn with 
the pressure dependence of Tn of pure CePdAl, we de- 
termine the effect of unit-cell-volume changes by exter- 
nal hydrostatic pressure p and by chemical pressure. 
We estimate the compressibility k = — {1/V) (dV/dp) 
of CePdAl from the bulk moduli By of the pure con- 
stituents (Voigt-Reuss-Hill approximation) as described 
in Ref. [SJI, yielding By = 70.8 GPa for CePdAl. Di- 
rect measurements of the compressibility of CePdAl per- 
formed at room temperature yield By = 47.8 GPa [35|. 
This value is rather low for HF systems. Furthermore 
the same authors also report on the pressure dependence 
of Tn where Tn is found to remain constant up to 0.3 
GPa and then decreases at a smaller rate than found in 
Ref. [l^l . We therefore prefer to employ the Voigt-Reuss- 
Hill approximation. We calculate the unit-cell volume 
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FIG. 3. Neel temperature Tn vs. unit-cell volume V of 
CePdi-xNixAl deduced from our data. The data for CeP- 
dAl under hydrostatic pressure were taken from Goto et al. 
f2^, (A: magnetic susceptibility, □: electrical resistivity, O: 
specific heat), see text for details. 



V(p) for each pressure by Goto et al. [22] to obtain 
Tn(V) (open symbols in Fig. [3]). The unit-cell volume 
V{x) of CePdi-icNiajAl was calculated from our x-ray 
diffraction data. The resulting Tn(V^) obtained from the 
data sets Tj^{x) and Tn(p) can be inferred from the lower 
horizontal scale of Fig. [3] We find that substituting Pd 
by Ni yields the same Tn (V) dependence as applying hy- 
drostatic pressure within the uncertainty of the measure- 
ments. This observation suggests that the decrease of Tn 
with Ni alloying or under pressure is a purely geometric, 
i.e., volume-dependent, effect. Thus, neither is disorder a 
driving factor for the suppression of magnetic order, nor 
can the decrease of Tn be attributed to a large quanti- 
tative change of the frustration in this system, since in 
these cases we would expect strongly different dependen- 
cies for Tn(p) and TN(a:). 

For the discussion of the possible effect of frustration 
on the QCP of CePdAl, the main problem is how to quan- 
tify the degree of frustration. Above we have already 
dismissed the Tn{V) dependence on the grounds that 
this dependence on the unit-cell volume is independent of 
whether V is changed by hydrostatic pressure or substi- 
tution of Pd by Ni. In geometrically frustrated insulating 
magnetic systems with stable moments, frustration arises 
from (near) cancelation of competing exchange interac- 
tions at the site of a given moment. In this case the frus- 
tration parameter ist defined a.s f — |6cw|/21si, where 
Gcw is the Curie- Weiss temperature derived from the 
susceptibility in the paramagnetic regime. Using this def- 
inition, it could be shown for the poorly metallic RXCui 
system with R = rare-earth element and A" = In or Cd 
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that there is an approximate inverse relation between / 
and the electrical conductivity a, i.e., higher a leads to 
lower / (Ref. [s^). If one were to take the same definition 
of / at face value for CePdi_a;Nia;Al, we would have to 
infer a strong increase of / with x because Tn decreases 
and |Ocw| remains roughly constant in the x range inves- 
tigated here, as is pictured in Fig. |4l Thus, the increase 
of the ratio |Ocw| /T^ with x can not be taken as indi- 
cator of increased geometric frustration, it rather reflects 
the competition between Kondo effect and RKKY inter- 
actions. A similar apparent increase of |6cw|/2n can 
be infered in CeCug-icAua;, which is a system without 
geometric frustration. 

Fig.|3]serves to illustrate an additional important find- 
ing. In the paramagnetic state the effective moment in- 
ferred from the slope of 1/x versus T for 30 < T < 300 K 
is practically the same for a; = and x = 0.144. The same 
holds for the intermediate concentrations [s^] . The inset 
of Fig. m reveals that, at low T, a marked difference in 
the T dependence develops. While x~^{T) varies super- 
linearly for x ~ Q above Tn, the data for x — 0.144 tend 
to show an exponent smaller than 1. This might possi- 
bly point to an anomalous scaling exponent a ~ 0.9 as 
observed for CeCue-iAua; [1, 0] although neutron scat- 
tering experiments are necessary to check this conjecture. 

What - then - might be the role of frustration concern- 
ing the QCP in CePdAl? We first discuss this issue in 
terms of the HMM model. While the magnetic structure 
of CePdAl is three-dimensional in nature with the order- 
ing wave vector Q = Or), our measurements of the 
specific heat C/T oc ln(To/T) for x ~ Xc and the lin- 
ear Tn(cc) dependence are indeed in accord with 2D AF 
quantum fluctuations in the HMM scenario. 

While inelastic neutron scattering experiments are nec- 
essary to identify possible 2D planes, our data can al- 
ready exclude certain possibilities. For a single kagome- 
like plane, the magnetic structure can be described by 
ferromagnetic chains, as mentioned above, that couple 
only weakly antiferromagnetically by the nun interaction 
J2, because of the frustration of the nn interactions Ji 



at the Ce(2) sites [26 



27l |. Taking into account the in- 
terlayer exchange J3, the magnetic order can be viewed 
as a stacking of ferromagnetic chains staggered along the 
c-direction interrupted by the frustrated Ce(2) moments 
thus leading to 2D ferromagnetic planes roughly perpen- 
dicular to the magnetic ordering wave vector Q. How- 
ever, 2D ferromagnetic fluctuations will give rise to a 
specific heat C/T cx T^a, which is incompatible with 
our data. A triple-Q structure is compatible with the 
neutron-scattering data if r = i. In this structure, 
however, AF planes cannot easily identified. In addi- 
tion the magnetic component t along the c direction 
is, in fact, incommensurate and varies slightly with T 
(t ^ 0.3535 -0.3585 between 0.35 and 2.7 K) While 
an incommensuration with small difference of r from i 
might in principle arise from stacking faults, the tem- 
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FIG. 4. Inverse susceptibility (x — Xo)~^ of CePdAl and 
CePdo.856Nio.i44Al versus temperature T. The inset shows 
the low temperature region below 30 K, where a distinct dif- 
ference in curvature emerges when approaching the QCP. A 
sample-dependent residual background contribution (see text) 
was subtracted from the data. The yellow line shows a fit with 
an power law l/x + 6 oc y"-^. 



perature dependence of r is not easy to explain within 
a triple-Q structure. Hence within the HMM model, 2D 
fluctuations conflned to AF planes perpendicular to the 
ah plane and separated by frustrated Ce(2) moments, are 
a viable explanation for the C/T cx log (Tq/T) behavior. 

While in the classical case of an anisotropic 3D system 
geometric frustration might give way to strictly 2D or- 
der, in a quantum system zero-point fluctuations might 
restore the original dimensionality [SS^. This, however, 
is not always the case. For example, it was shown in 
the context of a Bose-Einstein condensation of magnetic 
triplets in BaCuSi206 [39*1 that geometrical frustration 
yields a reduction of the spatial dimensionality [381, \^ . 
The partial magnetic frustration in CePdAl appears to 
reduce the dimensionality by forming frustrated planes 
from one third of the Ce moments, decoupling the an- 
tiferromagnetic planes of the other Ce moments. Thus 
CePdAl might be a unique system where frustration pro- 
vides a rationale how a system might find its way from 3D 
magnetic ordering to 2D criticality. Thus CePdAl might 
also be a very promising candidate to test the recent pre- 
dictions for the 2D-3D crossover in terms of the critical 
quasiparticle theory |4lJ . The deviations of Tn (x) from a 
linear x dependence close to Xc are compatible with such 
a crossover. 

Of course geometric frustration might give rise to ad- 
ditional fluctuations that are not included in the HMM 
model. Hence other scenarios might be conceivable to 
explain the quantum critical behavior of CePdAl. In- 
deed the susceptibility exponent a < 1 for a; = 0.144 
could possibly be due to anomalous scaling. Detailed 
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experiments on scaling relations between the energy, 
temperature and field dependencies as done, e.g., for 
CeCug-ajAua; 0, are necessary to narrow down the ori- 
gin of the unusual quantum critical behavior of CePdAl. 
As a further interesting point, CePdAl offers the possi- 
bility of investigating the Kondo effect of geometrically 
frustrated metals in detail. 

In conclusion, we have investigated systematically 
the effect of chemical pressure on the specific heat of 
CePdi_j;Nia;Al alloys for Ni concentrations ranging up 
to X — 0.144, where the magnetic order is completely 
suppressed, giving rise to a QCP. At the concentration- 
tuned QCP, the log (To /T) dependence of the specific 
heat suggests, if interpreted within the HMM scenario 
of conventional quantum criticality, the evolution of 2D 
magnetic fluctuations. The geometric frustration in this 
system leading to one third of the Ce moments not par- 
ticipating in long-range AF order, may thus be instru- 
mental in establishing the unusual character of quantum 
criticality, despite the fact that the magnetic order it- 
self is 3D. On the other hand, the partial order might 
lead to novel types of excitations close to the QCP not 
contained in the HMM model that cause the log T be- 
havior of C/T. Inelastic neutron scattering experiments 
are planned to test this conjecture. Further investiga- 
tions on single-crystalline samples of CePdi_a;Nij;Al at 
different magnetic fields applied along different crystal- 
lographic directions are under way in order to clarify the 
effect of anisotropy on quantum criticality in this par- 
tially frustrated system. 

We thank J. Schnialian and O. Stockert for fruitful 
discussions and S. Drobnik for help with some mea- 
surements. This work was supported by the Deutsche 
Forschungsgemeinschaft through FOR 960. H.v.L. 
thanks the National Science Foundation for partial sup- 
port of this work under Grant No. PHYS-1066293, and 
the Aspen Center for Physics for its hospitality. 



[1] 
[2] 
[3] 
[4] 

[5] 
[6] 

[7] 
[8] 



S. Sachdev, Quantum phase transitions (Wiley Online 
Library, 2007). 

S. Sondhi, S. Girvin, J. Carini, and D. Shahar, Reviews 
of Modern Physics 69, 315 (1997). 

H. V. Lohneysen, A. Rosch, M. Vojta, and P. Wolfle, Rev. 
Mod. Phys. 79, 1015 (2007). 

O. Stockert, H. v. Lohneysen, A. Rosch, N. Pyka, and 
M. Loewenhaupt, Phys. Rev. Lett. 80, 5627 (1998). 
A. Schroder, G. Aeppli, E. Bucher, R. Ramazashvili, and 
P. Coleman, Phys. Rev. Lett. 80, 5623 (1998). 
A. Schroder, G. Aeppli, R. Coldea, M. Adams, O. Stock- 
ert, H. Lohneysen, E. Bucher, R. Ramazashvili, and 
P. Coleman, Nature 407, 351 (2000). 
P. Coleman, C. Pepin, Q. Si, and R. Ramazashvili, J. 
Phys.: Cond. Mat. 13, R723 (2001). 
Q. Si, S. Rabello, K. Ingersent, and J. Smith, Nature 
413, 804 (2001). 



[9: 

[10 
[11 

[12; 
[13 



[14; 
[is: 

[16: 

[17 

[is: 

[19 



[20 
[21 

[22: 
[23: 

[24: 

[25: 

[26: 

[27 

[28: 

[29 

[3o: 

[31 

[32: 

[33 
[34: 

[35: 
[36 



J. A. Hertz, Phys. Rev. B 14, 1165 (1976). 
A. J. Millis, Phys. Rev. B 48, 7183 (1993). 
T. Moriya and T. Takimoto, J. Phys. Soc. Jpn. 64, 960 

(1995) . 

M. Garst, L. Fritz, A. Rosch, and M. Vojta, Phys. Rev. 
B 78, 235118 (2008). 

J. Custers, K. Lorenzer, M. Miiller, A. Prokofiev, 
A. Sidorenko, H. Winkler, A. Strydom, Y. Shimura, 
T. Sakakibara, R. Yu, et al., Nature Materials 11, 189 
(2012). 

P. A. Lee, Science 321, 1306 (2008). 

Q. Si, Physica B 378-80, 23 (2006). 

M. Vojta, Phys. Rev. B 78, 125109 (2008). 

P. Coleman and A. Nevidomskyy, J. Low Temp. Phys. 

161, 182 (2010). 

M. Kim and M. Aronson, J. Phys.: Condensed Matter 
23, 164204 (2011). 

A. Oyamada, K. Kamioka, K. Hashi, S. Maegawa, 
T. Goto, and H. Kitzawa, J. Phys. Soc. Jpn. 65, 128 

(1996) . 

A. Donni, G. Ehlers, H. Maletta, P. Fischer, H. Kitazawa, 
and M. Zolliker, J. Phys.: Cond. Mat. 8, 11213 (1996). 
Y. Isikawa, T. Mizushima, N. Fukushima, T. Kuwai, 
J. Sakurai, and H. Kitzawa, J. Phys. Soc. Jpn. 65 Suppl. 

B. 117 (1996). 

T. Goto, S. Hane, K. Umeo, T. Takabatake, and 

Y. Isikawa, J. Phys. Chem. Sol. 63, 1159 (2002). 

K. Prokes, P. Manuel, D. Adroja, H. Kitazawa, T. Goto, 

and Y. Isikawa, Physica B: Condensed Matter 385, 359 

(2006). 

A. Oyamada, S. Maegawa, M. Nishiyama, H. Kitazawa, 
and Y. Isikawa, Phys. Rev. B 77, 064432 (2008). 
P. Villars and K. Cenzual, Pearson's Crystal Data - Crys- 
tal Structure Database for Inorgamc Compounds ( on CD- 
Rom), Release 2011/12, ASM International, Materials 
Park, Ohio, USA. 

M. Niifiez-Regueiro, C. Lacroix, and B. Canals, Physica 
C: Superconductivity 282, 1885 (1997). 
Y. Motome, K. Nakamikawa, Y. Yamaji, and M. Uda- 
gawa, Phys. Rev. Lett. 105, 036403 (2010). 
Y. Isikawa, T. Kuwai, T. Mizushima, T. Abe, G. Naka- 
mura, and J. Sakurai, Physica B 281, 365 (2000). 
N. Bagrets, Ph.D. thesis, Universitat Karlsruhe (TH), 
Physikalisches Institut, Karlsruhe, Germany (2009). 
High-purity rare-earth metals acquired from Ma- 
terials Preparation Center, Ames Laboratory, US 
DOE Basic Energy Sciences, Ames, lA, USA, 
<|http://www.mpc.ameslab.gov >. 

A. Gribanov, A. Tursina, E. Murashova, Y. Seropegin, 
E. Bauer, H. Kaldarar, R. Lackner, H. Michor, E. Roya- 
nian, M. Reissner, et al., J. Phys.: Cond. Mat. 18, 9593 
(2006). 

C. Schank, F. Jahrling, L. Luo, A. Grauel, C. Wassilew, 
R. Borth, G. Olesch, C. Bredl, C. Geibel, and F. Steglich, 
J. All. & Comp. 207, 329 (1994). 

H. V. Lohneysen, J. Phys.: Condens. Matter 8, 9689 
(1996). 

G. Grimvall, Thermophysical properties of materials 
(The Royal Institute of Technology, Stockholm, Sweden, 
1999). 

J. Tang, A. Matsushita, H. Kitazawa, and T. Matsumoto, 
Physica B: Condensed Matter 217, 97 (1996). 
V. Fritsch, J. D. Thompson, J. L. Sarrao, H.-A. Krug von 
Nidda, R. M. Eremina, and A. Loidl, Phys. Rev. B 73, 



6 



094413 (2006). 

[37] V. Fritsch, C.-L. Huang, N. Bagrets, K. Grube, S. Schu- 
mann, and H. v. Lohneysen (2012), submitted to Phys. 
Stat. Sol. (b). 

[38] J. Schmalian and C. Batista, Phys. Rev. B 77, 094406 
(2008). 

[39] S. Sebastian, N. Harrison, C. Batista, L. Balicas, 



M. Jaime, P. Sharma, N. Kawashima, and I. Fisher, Na- 
ture 441, 617 (2006). 

[40] C. Batista, J. Schmalian, N. Kawashima, P. Sengupta, 
S. Sebastian, N. Harrison, M. Jaime, and I. Fisher, Phys. 
Rev. Lett. 98, 257201 (2007). 

[41] E. Abrahams and P. Wolfle, Proceedings of the National 
Academy of Sciences 109, 3238 (2012). 



